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ABSTRACT 

Although soft X-ray absorption features in AGN were discovered almost ten years ago, the nature 
and location of the gas creating them has remained controversial. Guided by the results of recent high- 
resolution X-ray spectroscopy, we argue that these features are created in a multi-temperature wind 
whose source of matter is photoionized evaporation from the inner edge of the obscuring torus often 
found surrounding AGN. Photoionized evaporation in the presence of a copious mass source locks the 
ratio of ionizing intensity to pressure to a critical value. However, a broad range of temperatures can all 
coexist in equilibrium for this value of the ratio of ionizing intensity to pressure. Consequently, the flow 
t-H \ is expected to be strongly inhomogeneous in temperature. The inferred distance of this material from the 

source of ionizing radiation depends on how much matter exists at the highest-obtainable temperature. 
This distance can be measured by monitoring how ionic column densities respond to changes in the 
^ , ionizing continuum on timescales of days to years. 

iy~j \ Subject headings: Galaxies: Active — Galaxies: Nuclei — Galaxies :Seyfert — Galaxies: Quasars: 

ff") . Absorption Lines — X-rays: Galaxies — Ultraviolet: Galaxies 

Q ' 1. INTRODUCTION 

X-ray spectra of active galactic nuclei (AGN) commonly reveal intrinsic absorption by highly ionized or "warm" gas. 
First seen in ROSAT and ASCA spectra (Turner et al. 1993; Reynolds et al. 1995; Reynolds 1997; George et al. 1998b), 
,-£< absorption edges of O VII and O VIII can be found in half or more of all type 1 Seyfert galaxies, and in a few quasars 
Dh| (Mathur et al. 1994). 

The location of this gas is still unknown. Suggestions range all the way from the broad- line region (BLR), possibly as 
;— i ' close as ~ 0.01 pc from the nucleus, to distances of 10 pc or more. Rapid apparent column density variations have been 
5/5 \ interpreted as suggesting locations in the BLR (Reynolds et al. 1995; George et al. 1998b). On the other hand, the same 
C$ data also indicate that at least some of the gas may be much farther away (Otani et al. 1996). Others have placed the gas 
at ~ 1 pc in the scattering region posited by Seyfert galaxy unification schemes (Krolik & Kriss 1995), or a factor of 10 
farther away (Bottorff et al. 2000). Morales et al. (2000) suggested that the absorber is spread all the way from < 0.05 pc 
to > 1 pc. Still less is known about the origin, dynamics, or destiny of the gas producing warm absorber features. Some 
of the diverse suggestions include: evaporation off "bloated stars" in the BLR (Netzer 1996); gas, evaporated off the torus 
obscuring the nucleus, that becomes the scattering gas seen in type 2 Seyfert galaxies (Krolik & Kriss 1995); and a wind 
driven off the accretion disk (Elvis 2000; Bottorff et al. 2000). 

With the launch of Chandra and XMM-Newton, we now have access to data of far better quality, and may reasonably 
hope to settle some of these outstanding questions. Most notably, the high resolution available with the grating spec- 
trographs on these spacecraft has permitted detection of resonance line features (Kaastra et al. 2000; Kaspi et al. 2000; 
Kaspi et al. 2001; Lee et al. 2001; Collinge et al. 2001). As we will show later, measurement (or even bounds) on line 
widths and shifts (i.e., relative to the host galaxy rest-frame) can provide powerful hints toward answering many of the 
questions just posed. It is our object here to show how this new data points us toward a revision of the "scattering 
wind" model. Because this model makes specific predictions about the character of variability on measurable timescales, 
it should be easily testable in the near future. 

2. THE MODEL: A CRITICAL-PRESSURE WIND 

2.1. Phenomenological guidance 

Measured column densities of different ionization stages seen in absorption can be used to constrain photoionization prop- 
erties of the absorbing gas under the assumption that the species involved are in ionization (and optionally also thermal) 
balance. The relative abundances of different species fix the ionization parameter £ = Li on /nr 2 ; the strengths of the lines 
and edges fix the hydrogen column density N . Here the ionizing luminosity Li on is conventionally taken to be the luminos- 
ity between 1 and 1000 Ryd, n is the hydrogen volume density, and r is the distance from the nucleus to the photoionized 
gas. Other, closely-related, versions of the ionization parameter are also used, e.g., U oxygen = L 53g de(L e /4irr 2 enc) 
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2 Warm Absorbers in AGN 

(George et al. 2000), where photon energy e is measured in keV and L e is the luminosity per unit energy. Its precise 
relation to £ depends on the shape of the continuum spectrum. 

These quantities can then be used (Turner ct al. 1993) to estimate the fractional radial thickness occupied by absorbing 
gas 

Ar/r = ZNr/L lon . (1) 

Placing the gas close to the nucleus implies a highly "clumped" structure, whereas larger distances imply a more volume- 
filling configuration. The gas cannot, in any case, be farther than r max = Li on /(N^) because Ar/r < 1. Simple 
photoionization-fitting cannot distinguish between a clumpy and a smooth configuration; it can only define the link 
between the degree of dumpiness (or smoothness) and the distance from the nucleus to the absorbing gas. 

The most commonly observed features (especially in ASCA data) are those associated with O VII and O VIII (Reynolds 1997; 
George et al. 1998b). In order to make them reasonably abundant, £ ~ 10 - 100. For this reason, many photoion- 
ization models of warm absorbers have suggested values of £ in this range, e.g., (Brandt et al. 1997; Reynolds 1997; 
George et al. 1998b; Mathur et al. 1997). On the other hand, more recent work has detected many other species of 
higher ionization level (Collinge et al. 2001), suggesting that there may be substantial amounts of gas with £ ~ 100 - 
1000 (Kaastra et al. 2000; Kaspi et al. 2001). The maximum distance of the absorbing region in a Seyfert galaxy is then 
Tmax ~ 30Lj on] 447V^ 1 ^{ Q pc, where the fiducial quantities (10 44 erg s _1 for Li on , 10 22 cm~ 2 for N, and 100 for £) have 
been chosen to be representative of those commonly inferred. 

Line profiles convey additional information. For example, Kaspi ct al. (2001) find that the absorption lines they detect 
have a mean blueshift relative to the galaxy's systemic velocity of 610 km s _1 . The gas must, therefore, possess a net 
flow toward us. The associated mass outflow rate is 

M out = AirCfinNrv = 0.72CN 22 r pc (v/500 km s" 1 ) M Q yr" 1 , (2) 

where its covering fraction around the nucleus is C, /in is the mean mass per H atom, r pc is the radius normalized to 
1 pc, and v is the mean outflow speed. Interestingly, the ratio of the outflow rate to the accretion rate required to power 
the nuclear luminosity can be ~ 100: 

Mout „ nnr , f11acc\ ( r \ ( L ion /L bol \ ( t, X~ ,„,, 

^- = 600C 5 oo (^-) [j- ) [—^- ) [ Wq ) , (3) 

where Ltoi is the bolometric luminosity and r) acc is the accretion efficiency (cf. Krolik & Begelman 1986, Reynolds 1997). 

The widths of the absorption lines appear to be similar to their blueshifts; many absorption lines also have emission 
components to the red (i.e., "P Cygni"-typc profiles) (Kaspi et al. 2001). The fact that the emission components are 
centered at the systemic velocity and have widths comparable to the blueshifts is most simply understood as indicating 
that emission and absorption are made by gas with similar properties, but that the line emission sources are located to 
the side of and behind the continuum source. 

All measures of gas velocity are roughly an order of magnitude smaller than the typical speeds of gas in the BLR. 
Following the common rule-of-thumb that speeds diminish at increasing distance from the central black hole, this contrast 
would suggest that the absorbing gas is rather more distant than the BLR (Kaspi et al. 2000). Still more speculatively, 
if velocities scale oc r -1 ' 2 , the warm absorber material must be ~ 100 times farther from the nucleus than the BLR, 
i.e., in NGC 3783 it must be parsecs from the nucleus if we adopt the reverberation-mapping estimate of its BLR scale 
(Reichcrt et al. 1994). 

On the basis of this evidence (modest velocity ouflow, maximum distance ~ 30Li on ^4 pc, moderately high ionization) 
and the comparisons to observations presented in §3, we find the most natural model to be a close relation to the one 
proposed by Krolik & Kriss (1995). Again, we identify the warm absorber with the hot scattering gas of Seyfert galaxy 
unification models. However, as discussed in the next subsection, our view of its thermodynamic state is somewhat 
different from the one expressed six years ago. 

2.2. Thermodynamics of the warm absorber 

In this model, the origin of the warm absorber gas is material evaporated off the inner edge of the toroidal ob- 
scuration. The position of this inner edge is determined by a combination of dust sublimation and photoionization 

(Pier & Voit 1995), but is expected to be ~ ^L^ n 44 pc. The existence of rings of H2O masers at approximately this 
distance (Greenhill et al. 1996) is in keeping with this expectation. Thus, the source for the warm absorbing gas would 
be roughly this far from the nucleus. 

The thermodynamics of this evaporating gas are best described by an ionization parameter framed in terms of the gas 
pressure, such as S = L ion /(4Trr 2 cnkT) (Krolik, McKee & Tarter 1981). When S exceeds a critical value (S c ^10-30 
for typical AGN spectra), the gas is no longer capable of maintaining a cool (~ 10 4 K) equilibrium temperature; instead 
its equilibrium temperature rises very rapidly with increasing S until it approaches the Compton temperature, which in 
AGN is generally several orders of magnitude greater. 

The detailed character of this rise depends on the specific shape of the continuum (and also on details of atomic physics 
that are imperfectly known), but its qualitative behavior — a nearly vertical rise from temperatures ~3x 10 4 K to ~ 10 6 K 
at S ~ S c — is the same for almost any AGN-like spectrum. In Figure 1 we show T(E) (as computed using Version 2.1 of 
XSTAR, Kallman 2000) for two continuum shapes based on data from the type 1 Seyfert galaxy NGC 3783, one suggested 
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Fig. 1. — The equilibrium temperature as a function of 
S for the Kaiser et al. (2001) spectrum of NGC 3783 
(solid curve) and the Kaspi et al. (2001) spectrum 
(dashed curve). 
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Fig. 2 — Proposed continuum shapes for NGC 3783: 
solid curve by Kaiser et al. (2001), dashed curve by 
Kaspi et al. (2001). Both are plotted in terms of flux 
per logarithmic frequency interval, vF v ] the vertical 
scale is arbitrary. 



by Kaiser et al. (2001), the other by Kaspi et al. (2001). The rise is very similar for both spectra; almost the only 
difference is the offset in S c . The fiducial £ = 100 corresponds to a temperature part way up this rise, 1.5 x 10 5 K for the 
Kaiser et al. spectrum, 5.5 x 10 4 K for the Kaspi et al. spectrum. In terms of £, the vertical portion of the equilibrium 
curve runs from £ ~ 10 to £ ~ 1000 in the case of the Kaiser et al. spectrum, but from £ ~ 40 to £ ~ 1500 for the Kaspi 
et al. spectrum. 

In this paper we are primarily interested in material near and outside the obscuring torus. In keeping with theoretical 
expectations that the torus should reradiate a substantial portion of the nuclear luminosity in the infrared, we have 
required that in each spectrum the luminosity at energies below 2 eV match the ionizing luminosity. No alteration in the 
Kaiser et al. spectrum was necessary to accomplish this, but to achieve this with the Kaspi et al. spectrum we added a 
new component at 0.1 eV. Both spectra are shown in Figure 2. The additional infrared intensity makes little difference 
to the equilibrium curve at temperatures below ~ 10 6 K; the primary effect is to reduce the Compton temperature and 
create the "jog" toward higher 5 in the range 1 - 3 x 10 6 K. Closer to the nucleus, the infrared intensity would be relatively 
weaker, leading to a higher Compton temperature and a longer section of nearly-vertical temperature rise. Throughout 
the remainder of this paper, specific numbers derived from photoionization calculations assume the Kaiser et al. spectral 
shape. 

Places where the slope dT/dE > are stable against isobaric thermal perturbations, but places where the curve "doubles 
back" (i.e., where dT/dE < 0) are thermally unstable. As shown by the curves of Figure 1, near the obscuring torus there 
is a significant stretch of the equilibrium curve where the sign of dT/dE is ill-defined because the temperature rises almost 
exactly vertically: dE/dT ~ for 10 45 K^ T <; 10 59 K (and \dE/dT\ remains very small up to T ~ 10 65 K). This region 
is therefore better described as having approximately marginal stability. That is, thermal fluctuations neither grow nor 
decay exponentially; there is little feedback in either sense. Although the equilibrium curves are never exactly vertical 
locally, because the overall sense is very nearly vertical, tiny pressure perturbations in regions that are nominally stable 
can nonetheless lead to large temperature changes. Likewise, temperature perturbations in regions of nominal instability 
do not lead to arbitrarily large temperature changes — at constant pressure, material rises or falls along a line of constant 
E, and cannot go far before intersecting the equilibrium curve again, this time on a locally stable segment. It is in this 
sense that we describe this portion of the equilibrium curve as being "marginally stable" . 

The reason to pay such close attention to the character of equilibria with E ~ E c is that when there is a copious source 
of gas available for rapid heating (as the obscuring torus provides in this model), rather than being a special value, E c is 
a preferred value of the ionization parameter. Suppose that the pressure in the gas at the illuminated face of the torus 
were very low, so that S > E c . Then this gas would be subject to strong net heating. As its temperature rises, it expands 
and its opacity drops sharply. Deep inside the torus, the large continuum optical depth protects the gas and keeps it 
relatively cool (the most likely temperature is <J 10 3 K: Krolik & Lepp 1989, Neufeld et al. 1994), but as the opacity of 
matter at the illuminated edge is eliminated, new gas is exposed and becomes available for "evaporation" . As more and 
more gas is removed from the torus and added to the hot phase, the pressure of the hot phase rises because (as we shall 
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show shortly), its heating time is much shorter than the flow time. Eventually, the pressure rises high enough that S 
falls to S c . Conversely, if the ambient pressure were so high that S < 5 C , gas outside the torus would cool and condense, 
diminishing the pressure. 

Thus, where the gas is injected, the ionization parameter is fixed at 5 C . Moreover, with one relatively modest proviso, it 
remains at that value even as the gas flows away from its injection point. So long as its ionization and thermal equilibration 
times arc short compared to the flow time, the gas must follow the equilibrium curve. Because the equilibrium curve lies 
along the vertical line S = S c until the temperature exceeds ~ 10 6 K, the gas's ionization parameter remains fixed at 
that value until the majority of its mass reaches high temperatures. 

In fact, both the ionization and thermal balance timescales are indeed considerably shorter than the flow time over 
a wide range of conditions. The flow timescale tfi ow = r/v ~6x 10 10 r P c(f/500 km s -1 ) -1 s. On the other hand, the 
thermal equilibration time t coo i ~2x 10 7 s for our fiducial parameters (£ = 100, T = 1.5 x 10 5 K) when the gas is placed 
at 1 pc. In Table 1 we list abundances and ionization equilibration timescales ti on for commonly observed ions, calculated 
at ionization parameters £ = 30, 100, 300, and 1000 for gas illuminated at 1 pc by a central source with an ionizing 
luminosity of 10 44 erg s _1 . For these parameters, the most abundant ions have t, „ ~ 10 6±0 - 5 s. Roughly speaking, 
tcool, ion ft flow f* r , so thermal and ionization equilibrium both become better approximations when the gas is clumpier 
and closer to the nucleus, but even at r max thermal and ionization equilibrium should prevail for a substantial range of 
temperatures at S = S c . 

At fixed radius and pressure, both t coo i and ti on vary along the vertical portion of the T"(S) curve. Roughly speaking, 
the thermal timescale grows ex T 2 at constant S. At the highest temperatures, however, t coo i grows more rapidly as line 
cooling is replaced by Compton cooling. The Compton equilibration time is long enough, ~ 10 12 r 2 c L^l s, that fully 
reaching Comptonization equilibrium is unlikely 1 . Moreover, when t coo i ~ tfi ow , adiabatic cooling can become important 
(as in Krolik & Begelman 1986 and Balsara & Krolik 1993). Thus, the highest temperature on the equilibrium curve 

reachable by gas injected cool is the temperature at which t coo i < tfi ow , here T max ~ 4 x 10 6 rp C ' K (for distances J> 1 pc; 
at smaller distances, T max approaches the Compton temperature). When the gas temperature is <J T max , it is no longer 
constrained to lie on the equilibrium curve; due to adiabatic cooling, it tends to move into the region of net heating. 

The ionization equilibration timescale (at fixed S but variable £ or T) behaves in a somewhat more complicated fashion. 
It is determined by a combination of the ionizing flux and the local density. The former dominates when the ionization 
rate for ion j exceeds the recombination rate from j to j — 1, the latter when the sense of this comparison reverses. When 
the most abundant ion of a given element has charge > j, the former case generally applies. Because many of the species 
most important for soft X-ray features arc H- and He-like and these elements are often mostly stripped, the former case 
applies more often than the latter. 

Although the flux at fixed photon energy is, of course, entirely independent of T at a given radius, as the gas changes 
temperature, the photon energies relevant to its ionization state change. Both the photon fluxes and the cross sections 
can be strong functions of energy. For this reason, ti on also generally rises with increasing temperature, but rather more 
slowly than t coo i — for the dominant unstripped ion of each element, ti on increases by factors of ~ 1 - 5 from the bottom 
to the top of the vertical branch of the equilibrium curve. 

We conclude, then, that over a wide range of possible radii, the thermal and ionization equilibration times are indeed 
short compared to the flow time over most of the vertical section of the thermal equilibrium curve. After a flow time, it is 
possible for some gas to reach high enough temperatures that it is no longer required to be in thermal balance; however, 
it is not necessary that all the gas reach this state. Marginal thermal stability means that the temperature of a given gas 
parcel can easily move either up or down on timescales ~ t coo i; all that is required is that its density change enough to 
compensate for the temperature change. Thus, we expect a range of temperatures to be present in this region, and also 
a corresponding range of densities. 

2.3. Dynamics of the warm absorber 

These considerations also constrain the absorber's dynamics. The absorber's equation of state is determined by how 
much of its material is at low enough temperature to demand thermal balance. When most of the gas is in such a state, 
S is fixed at S c , and the pressure must vary oc r~ 2 . The strength of the pressure gradient relative to gravity is then given 

Here r\ ou t is the ratio of the mass loss rate to L/c 2 , and we suppose that r is comparable to the inner edge of the torus, 
where the radiation effective temperature is T c . The fiducial parameters in Equation 4 are 1O 6 M0 for the central mass 
M and 1000 K for T c (the latter chosen with reference to the results of Pier & Voit 1995). 

For there to be outward acceleration, Qu > (1 — kL/ Le), where u is the speed in free-fall units and k is the frequency- 
averaged opacity in Thomson units. This condition may be regarded as an upper bound on the mass loss rate in the sense 
that rj out <J 240u(L/Le) 1 ' 4 M 6 (assuming that radiation forces do not overcome gravity). Strikingly, this bound roughly 
matches the empirical estimate derived from eq. 3. In terms of the notation of that equation, rj out = M out / (M acc r) acc ) . 
Note, however, when making this comparision, that formally M out as estimated in equation 2 reflects only the mass loss 

1 Note that the Compton time depends on the total luminosity, not just the ionizing luminosity. 
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associated with matter at a single value of £; if there is a distribution (as we have now argued there should be), M out is 
the sum of the contributions from all the components (see also §3.2). 

The observed spectra suggest that the radiation pressure contribution is probably secondary. The ratio between the 
amount of photon momentum scattered by the wind and the total wind momentum is 

AL £r/Ar _ c in _ 3 ( AL/L\( £\ f _ , A ^ ( v \- 2 (N(£ = 100) 

100 



Jrad 



»*^m 



(r/Ar) 



(5) 



L Att^hcv 2 V 0- 1 / V 100 / 'VSOOkms- 1 / V N tot 

where the last ratio accounts for the range of ionization parameters that may be present. A scattered fraction AL/ L ~ 0.1 
appears to be consistent with what is seen (Kaspi ct al. 2001). Unless Ar/r is quite small, f ra d is likely to be well less 
than unity. 

The pressure in the wind is ex (T)M(u/r 1 * 2 )~ 1 r~ 2 , where (T) is the mean temperature. Fixed S then requires 
(T)M(u/r ' ) to be constant. If M were constant (which is not strictly required due to possible continuing evaporation 
or condensation), the velocity (and (T)) would increase only logarithmically with r. However, because the gas dynamics 
depend on its pressure, not its temperature, there is no requirement for the temperature to rise smoothly and homo- 
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X 
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2.45 x 


IO 5 
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x IO- 6 


2.46 x 


IO 5 


2.49 x IO- 7 


2.49 x IO 5 


N V 


9.76 x 


io- 7 


1.75 x 10 4 
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X 


io- 9 


1.75 x 


IO 4 
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1.66 x 


IO 4 


8.58 x IO- 14 


1.65 x IO 4 


N VI 


2.20 x 


io- 5 


1.30 x 10 5 


9.63 


X 


io- 7 


1.30 x 
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X 


io- 4 


4.21 x 


IO 5 
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4.23 x 


IO 5 
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io- 7 


3.04 x 


IO 4 
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x IO" 9 


2.98 x 
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10" s 
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X 
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6.80 x 
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io- 6 
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io- 5 


5.28 x 
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IO 5 
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NeX 


3.46 x 


io- 6 


1.49 x IO 6 
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X 


io- 5 


1.49 x 


IO 6 
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x IO- 6 


1.50 x 


IO 6 


2.01 x IO- 6 


1.52 x IO 6 


MgXI 


6.33 x 


io- 6 


1.08 x IO 6 
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X 


io- 5 


1.08 x 


IO 6 


3.43 


x IO- 6 


1.08 x 


IO 6 


2.21 x IO- 7 
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MgXII 


3.40 x 


io- 7 
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X 


io- 6 


2.85 x 


IO 6 
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x IO- 5 


2.86 x 


IO 6 


4.91 x IO- 6 


2.89 x IO 6 


Si XIII 


9.34 x 


io- 7 


1.98 x IO 6 


1.01 


X 


io- 5 


1.99 x 


10" 


9.50 


x IO- 6 


1.99 x 


IO 6 


1.09 x IO- 6 


2.02 x IO 6 


Si XIV 


1.27 x 


io- 8 


4.91 x IO 6 


1.03 


X 


io- 6 


4.90 x 


10" 


1.04 


x IO- 5 


4.92 x 


IO 6 


5.89 x IO- 6 


4.99 x IO 6 


sxv 


4.75 x 


io- 8 


3.29 x IO 6 


1.72 


X 


io- 6 


3.29 x 


IO 6 


9.44 


x IO- 6 


3.31 x 


IO 6 


2.22 x IO- 6 


3.35 x IO 6 


SXVI 


6.38 x 


L0- 11 


2.97 x 10 12 


1.82 


X 


IO -8 


7.89 x 


IO 6 


1.17 


x IO- 6 


7.92 x 


IO 6 


1.36 x IO- 6 


8.03 x IO 6 


Ar XVII 


3.32 x 


L0- 10 


5.03 x IO 6 


8.07 


X 


io- 8 


5.04 x 


IO 6 


1.51 


x IO- 6 


5.05 x 


IO 6 


7.71 x IO- 7 


5.12 x IO 6 


Ar XVIII 


9.05 x 


L0~ 13 


1.05 x IO 6 


1.85 


X 


io- 9 


1.20 x 


IO 7 


4.32 


x IO- 7 


1.20 x 


IO 7 


1.08 x IO- 6 


1.22 x IO 7 


CaXIX 


2.79 x 


L0- 12 


7.58 x IO 6 


3.85 


X 


io- 9 


7.58 x 


IO 6 


4.06 


x IO" 7 


7.61 x 


IO 6 


9.02 x IO- 7 


7.71 x IO 6 


CaXX 


8.36 x 


L0- 14 


2.32 x IO 3 


9.17 


X 


io- 10 


1.89 x 


IO 7 


1.19 


x IO- 6 


1.75 x 


IO 7 


1.29 x IO- 5 


1.77 x IO 7 


FeXVI 


6.35 x 


io- 6 


1.25 x 10 5 


1.43 


X 


io- 6 


1.25 x 


IO 5 


1.45 


x IO" 8 


1.26 x 


IO 5 


4.06 x IO- 13 


1.27 x IO 5 


Fe XVII 


4.35 x 


io- 6 


1.78 x 10 5 


9.40 


X 


io- 6 


1.78 x 


IO 5 


6.48 


x IO- 7 


1.79 x 


IO 5 


7.61 x IO- 11 


1.81 x IO 5 


Fe XVIII 


6.71 x 


io- 7 


2.57 x 10 5 


9.24 


X 


io- 6 


2.58 x 


IO 5 


4.33 


x IO- 6 


2.59 x 


IO 5 


3.73 x IO" 9 


2.62 x IO 5 


FeXIX 


5.66 x 


io- 8 


2.72 x 10 5 


3.81 


X 


io- 6 


2.72 x 


IO 5 


8.05 


x IO- 6 


2.73 x 


IO 5 


4.84 x IO -8 


2.76 x IO 5 


FeXX 


5.27 x 


IO- 9 


4.92 x IO 5 


1.32 


X 


io- 6 


4.93 x 


10 s 


1.08 


x IO- 5 


4.95 x 


IO 5 


5.50 x IO- 7 


5.01 x IO 5 


FeXXI 


1.26 x 


L0- 10 


3.86 x IO 5 


1.06 


X 


io- 7 


3.86 x 


10 s 


3.66 


x IO- 6 


3.88 x 


IO 5 


1.13 x IO- 6 


3.92 x IO 5 


Fc XXII 








4.05 


X 


IO -8 


9.96 x 


IO 5 


3.04 


x IO- 6 


1.00 x 


IO 6 


4.83 x IO- 6 


1.01 x IO 6 


Fc XXIII 








4.62 


X 


io- 9 


2.15 x 


IO 6 


8.39 


x IO -7 


2.16 x 


IO 6 


6.88 x IO- 6 


2.18 x IO 6 


Fe XXIV 








1.08 


X 


io- 10 


3.50 x 


10" 


1.00 


x IO- 7 


3.51 x 


IO 6 


5.09 x IO- 6 


3.55 x IO 6 


FeXXV 


















6.15 


x IO -8 


2.08 x 


IO 7 


1.07 x IO- 5 


2.11 x IO 7 


Fc XXVI 


















2.50 


x IO -9 


4.68 x 


IO 7 


2.20 x 1Q- 6 


4.64 x IO 7 



Note. — Abundances relative to hydrogen (A) and equilibration timescales (r oq ) for selected ions of the most astrophysically- 
abundant elements in thermal balance at photoionization parameters of £ = 30, 100, 300, 1000 assuming a distance of r = 1 pc from 
a source of luminosity IO 44 erg s _1 . The computations were done using XSTAR Version 2.1. 
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geneously so long as S = 5 C . Temperature fluctuations with wavelengths shorter than r and long enough that thermal 
conduction does not damp them can be easily created. We would expect, therefore, that, although there may be an overall 
sense of temperature evolution, there can be substantial local inhomogcneity. Eventually, the approximation made at the 
beginning that S = S c must break down. As Compton equilibrium is ncared, the equilibrium 5 begins to increase and 
the requirement of tracking the equilibrium curve is weakened as t coo i approaches and then exceeds tfi ow ■ 

3. COMPARISON WITH OBSERVATIONS 

3.1. Inferred ionization parameters and temperatures 

As we have just seen, this model predicts that the distribution of S in this gas should be very narrow, but the temperature 
distribution is likely to be much wider. In a fashion strikingly consistent with these predictions, Kaspi et al. (2001) found 
it impossible to adequately fit their data with only a single component. Moreover, their two inferred components lie 
squarely on the vertical section of the equilibrium curve, one near 8 x 10 4 K, the other at ~ 6 x 10 5 K. The similarity of all 
the line profiles suggests that these two components are subject to more or less the same forces, and are therefore not very 
distant from each other. Because they have the same value of S, if they lie at similar distances from the central source, 
their pressures must also be the same. Only if there is global pressure-confinement (i.e., a volume-filling background 
medium) would one expect these pressures to match. 

Our model predicts temperatures across a somewhat wider range, from ~ 3 x 10 4 K up to ~ 10 6 K. It may be difficult 
to detect the hottest gas, however. As one can see in Table 1, only for elements with atomic number > 12 (Mg and above) 
do more than 10% of the atoms retain at least one electron when £ = 1000 (T~9x 10 5 K); only a small increase in 
temperature beyond that point suffices to eliminate almost all the unstripped ions except Fe. In the band covered by the 
Chandra gratings, only two lines in gas with £ ~ 1000 (Fe XXII 10.98,11.77) can be expected to have line-center optical 
depths greater than unity when the velocity span is 500 km s^ 1 and the column density is 10 22 cm~ 2 (although another 
half dozen lines have opacities a factor of two smaller). Moreover, because the H-like stages of Mg, Si, S, Ar begin to be 
abundant when T>5x 10 5 K, it can be hard to distinguish 10 6 K gas from gas at half that temperature. 

UV features, on the other hand, grow in opacity at the lowest temperatures. For example, Table 1 shows that the 
fractional abundances 2 of C IV and O VI relative to H rise from 5.8 x 10~ 10 and 5.4 x 10 -7 , respectively, at £ = 100 
(1.5 x 10 5 K) to 2.4 x 10~ 7 and 3.9 x 10 -5 at £ = 30 (4.7 x 10 4 K). Any significant amount of gas at the low end of 
the temperature range is then likely to generate UV absorption features. These arc, in fact, often observed. In order to 
produce the line strengths seen, not much is required: for typical cases (Crenshaw et al. 1999; Kriss et al. 2000; Kaiser 
ct al. 2001), <; 10 19 H cm~ 2 suffices if the gas temperature is 5 x 10 4 K. 

3.2. Clumping 

If there is nothing to restrain it, there is no reason why this gas should not expand to fill the volume available to it. In 
fact, if the absorber gas were not volume- filling, one would have to ask what gas occupied the remainder of the region, and 
that gas would be subject to the same thermodynamic constraints as the gas under consideration. Within this volume, 
however, there are no obvious constraints favoring any particular temperature from ~ 3 x 10 4 K up to the maximum 
temperature achievable within a flow time, a few xlO 6 K. 

We do not know how to predict the distribution of column density with temperature dN/dT. Once it is given, though, we 
can find the distribution of volume with temperature: A differential column density dN having volume density n occupies 
a radial shell with thickness dAr = dN/n. In terms of the pressure p, dAr = (kT/p)dN, so that the distribution of 
fractional volume with temperature d(Ar/r)/dT oc TdN/dT. The normalization of this relation is fixed by the constraint 
J dT d(Ar/r)/dT = 1. That is, for equal column density and fixed pressure, higher temperature regions occupy larger 
volumes. 

Given a particular distribution dN/dT, it is possible to predict the distance of the entire system from the central 
continuum source. Suppose, for example, that dN/d\ogT is constant from T m i n = 3 x 10 4 K up to T max — 10 6 K. 
Because £ oc S/T, from Equation 1 it then follows that the distance at which this material exists is 

l-'icm *-Og\l max / -L min ) *j I ^tot \ / ^c \ / ^-raax \ /-.\ 

rpred ~ AircZ c kT max N tot ~ 4 ^°"< 44 v 3.5 x 10 22 cm" 2 ) [lb J {TVWkJ P °' W 

where N tot is the total column density, approximately a factor of 3.5 greater than the column density per logT. Unfortu- 
nately, this estimate is rather model-dependent because we have no way of determining the true dN/dT. In particular, to 
the degree that the higher temperatures have larger column densities, r prec i would decrease. However, this example does 
emphasize the likelihood that gas at intermediate T occupies only a fraction of the volume. 

3.3. Variability 

For considerations of bulk dynamics, ionization balance is a good approximation when ti on < tfi ow . However, because 
we can observe changes on timescales much shorter than tfiow, it is possible for us to see continuum fluctuations driving 
short-lived changes in the ionization balance. In fact, because the observed continuum flux and the absorption features 
"ride the light-cone" together, we can even follow fluctuations on timescales much shorter than the light-crossing time. 

2 XSTAR uses the elemental abundances of Grcvcssc ct al. 1996. 
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The column densities we measure reflect the history of the continuum flux over the past ti on ] the temperature inferred 
from column densities is the product of heating and cooling over the past t coo i . 

The relation between observed absorption column densities and the continuum provides an excellent diagnostic of the 
actual distance between the continuum source and the gas. At a distance of 1 pc, Table 1 shows ti on ~ 10 6±0 - 5 s for those 
species most visible in the soft X-ray band. However, since ti on oc r 2 , the magnitude of ti on depends sensitively on the 
column density in the hottest (and nearly transparent) component. For example, if the guess illustrated in equation 6 is 
correct, the true distance would be ~ 4 times larger, making the ionization timescales ~ 16 times longer, or l0 7 - 2±0 - 5 s for 
the species creating soft X-ray features. UV features should respond somewhat more quickly. Although it may be very 
difficult to monitor AGN X-ray lightcurves with sampling density sufficient to quantitatively match continuum fluctuations 
with ionization fluctuations, it should at least be possible to estimate upper bounds on ionization response times. These 
bounds can then provide constraints on the distance to this material, and, in this model, dN/dT. 

As of now, there has been little in the way of systematic monitoring of warm absorbers. The clearest report in the 
literature of very rapid changes in an absorbing column is an apparent sharp increase in the O VIII K-edge optical depth 
in MCG-6-30-15 that occurred in a span of ~ 10 4 s during an ASCA observation (Reynolds ct al. 1995; Otani et al. 1996). 
Within the context of our model, this time could be brought into agreement with the observed time if there were a column 
density ~2x 10 23 H cm~ 2 in gas at T ;> 10 6 K. However, it is also possible that this feature may have been confused with 
O VIII Ka emission broadened by relativistic motions in the accretion disk (Branduardi-Raymont et al. 2001). The fact 
that when the O VIII column density in MCG -6-30-15 increased, the measured O VII column density was constant may 
support this alternate interpretation. A somewhat similar event occurred in NGC 4051, but with the roles of O VII and 
O VIII reversed — the O VII column density changed while the O VIII column density was constant (Guainazzi ct al. 1996). 
In this case, because the luminosity of the nucleus is so small (~ 10 42 erg s _1 ), the timescale — ~ 10 4 s — places only a 
weak upper bound on the distance, r <J 0.2 pc. This is, of course, easily consistent with estimates like equation 6. 

Evidence in the UV to date is similarly fragmentary. Kriss et al. (1997) reported Lyman edge optical depth changes 
in NGC 4151, and Espey et al. (1998) saw C III column density changes in the same object that anti-correlated with the 
continuum flux on timescales of days. Unfortunately, the X-ray absorption in NGC 4151 is so complex that the warm 
absorber column density cannot be estimated well enough for us to predict the ionization timescale. 

3.4. Emission lines and the covering factor 

The peculiar geometry of flow through the opening in toroidal obscuration means that we cannot see the entire emission 
portion of P Cygni profiles created by absorption lines. In part this is because at many radii the warm ionized gas simply 
does not fill all solid angle (most of it is occupied by the dusty, molecular obscuring gas); in part this is because the 
obscuring matter partially blocks our view of the far side of the flow. In line with this prediction, Kaspi et al. (2001) use 
the P Cygni features they discovered in the spectrum of NGC 3783 to estimate covering fractions of ~ 0.3 and ~ 0.5 for 
their high- and low-ionization components, respectively. 

Unless r <C r max , variations in the column densities of the emission components on the timescale of continuum fluctua- 
tions (days to weeks) should be averaged out because the light-crossing time is much longer (several years) . Consequently, 
the covering fraction inferred on the basis of a single "snap-shot" may not be a fair measure of the true covering fraction. 

4. SUMMARY 

A wide variety of evidence now points toward locating the gas that produces ionized absorption features in the soft 
X-ray spectra of AGN in a warm wind that lies a few parsecs from the nucleus in Scyfert galaxies. Both absorption and 
emission lines are seen with widths of ~ 500 km s _1 , in many cases with the telltale "P Cygni" velocity shifts. The 
inferred ionization parameters and column densities indicate maximum distances ~ 30 pc. 

The specific ionization states indicated by the observed lines suggest that there are a number of components with 
different values of £, the density form of the ionization parameter. However, these all have the same value (S = 5 C ) of 
the pressure form of the ionization parameter. Above S c , low-temperature equilibria do not exist, and the temperature 
rises from ~ 3 x 10 4 K to ~ 10 6 K with S very nearly constant at E c . If these components coexist at the same distance 
from the central source, they share the same pressure, suggesting a volume-filling configuration by which the pressure is 
regulated. If the volume-filling phase is the highest attainable temperature on the marginally stable branch of the thermal 
equilibrium curve, the actual distance could be substantially smaller than the maximum; a factor of ten smaller would 
not be surprising. 

Because the thermal equilibria at 5 = S c are marginally stable to fixed-pressure thermal perturbations, material can 
move easily from one temperature to another. The only constraint is that accessible temperatures must have thermal 
equilibration times short compared to the flow time. Consequently, we expect that multi-temperature models will in 
general be required in order to explain all the observed features. 

Given that a very similar wind in a very similar location has been previously inferred from the scattering of nuclear 
light in obscured Seyfert galaxies, it makes sense to suppose that the two are essentially the same structure. However, 
absorption features in type 1 Seyfert galaxy X-ray spectra provide much more sensitive diagnostics of its state than does 
the polarized reflection of the continuum that can be observed in type 2 Seyfert galaxies. In particular, the time-dependent 
response of the absorbing column densities as the X-ray continuum level changes should provide a very sensitive test of 
these ideas. 
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